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I.  INTRODUCTION 

An  approximate  asymptotic  high  frequency  solution  is  developed  for  the 
field  exterior  to  a  smooth,  perfectly  conducting  convex  surface  when  it  is 
excited  by  a  ray  optical  electromagnetic  field.  In  the  conventional  GTD 
solution  to  this  problem,  the  total  field  exterior  to  a  surface  is  associated 
with  the  usual  geometrical  optics  (GO)  incident  and  reflected  rays  in  the 
lit  (or  illuminated)  region;  whereas,  in  the  shadow  region  it  is  associated 
entirely  with  the  surface  diffracted  rays*  introduced  by  Keller  [1,?,3]. 

These  rays  are  illustrated  in  Figure  1  together  with  the  regions  labeled 
I,  II,  III  and  V  which  divide  the  space  exterior  to  the  convex  surface  as 
follows.  The  shaded  region  II  constitutes  the  transition  region  adjacent 
to  the  shadow  boundary  (SB  of  Figure  1)  which  divides  the  space  exterior 
to  the  scatterer  into  the  lit  and  shadow  zones.  The  angular  extent  of 
region  II  is  of  the  order  of  (2/kPg(Q^))  '  radians,  where  k  is  the  wave- 
number  of  the  homogeneous,  isotropic  medium  exterior  to  the  convex  surface, 
and  Pg(Q])  is  the  radius  of  curvature  of  the  convex  surface,  in  the  plane 
containing  the  incident  ray  and  the  surface  normal.  The  point  of  grazing 
incidence  is  specified  by  Q^.  The  sub  regions  IV  and  VI  in  the  shadow  and 
lit  zones,  respectively,  denote  those  portions  of  region  II  which  are  ex¬ 
tremely  close  to  the  surface.  In  particular,  regions  IV  and  V  are  close 
to  the  portion  of  the  surface  which  is  a  caustic  of  the  surface  diffracted 
rays;  whereas,  region  VI  is  in  the  vicinity  of  which  is  a  caustic  of  the 
reflected  ray  for  grazing  incidence.  Thus,  regions  IV,  V  and  VI  are  commonly 
referred  to  as  the  caustic  or  surface  boundary  layer  regions.  The  GTD  is 
valid  in  regions  I  and  III,  but  it  fails  in  regions  II  and  V. 

*The  surface  diffracted  rays  are  in  general  present  even  in  the  illuminated 


region  if  the  surface  is  closed. 
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Pl  =  FIELD  POINT  IN  LIT  ZONE 
Pg*  FIELD  POINT  IN  SHADOW  ZONE 


Figure  1.  The  rays  and  regions  associated  with  scattering  by 
a  smooth  convex  cylinder. 


An  asymptotic  analysis  for  the  fields  in  regions  II  and  V  is  compli¬ 
cated  due  to  the  fact  that  the  asymptotic  character  of  the  field  changes 
continuously  but  rapidly  from  one  form  to  another  across  each  of  these  re¬ 
gions.  In  this  paper,  an  approximate  uniform  GTD  (orUGTD)  solution  is  de¬ 
veloped  which  remains  valid  in  regions  I,  III  and  the  portion  of  region  II 
which  lies  outside  the  sub  regions  IV  and  VI.  An  approximate  asymptotic 
solution  which  is  valid  within  the  caustic  boundary  layer  regions  IV,  V  and 
VI  is  not  considered  in  this  work;  it  will  be  reported  separately  in  the 
future. 

The  present  ansatz  or  formulation  of  the  uniform  GTD  solution  is  based 
on  an  asymptotic  solution  recently  obtained  by  PathakCA]  for  the  canonical 
problem  of  plane  wave  scattering  by  a  perfectly  conducting  circular  cylinder. 
The  work  in [4]  extends  and  improves  some  of  our  earlier  work  [5],  and  like 
most  other  analyses  dealing  with  the  diffraction  by  convex  surfaces,  this 
development  draws  upon  the  pioneering  contributions  of  Fock[6].  The  sim¬ 
ilarities  and  differences  between  the  solution  of  the  canonical  problem  in 
[4]  and  some  of  the  other  solutions  which  currently  exist  in  the  literature 
[7,8,9,10,11]  are  also  discussed  in  [4].  The  present  ansatz  based  on  the 
form  of  the  solution  of  the  canonical  problem  in  [4]  leads  to  a  uniform  GTD 
solution  for  an  arbitrary  convex  surface  which  is  convenient  for  engineering 
applications. 

It  is  assumed  here  that  the  incident  electromagnetic  fields  may  be 
approximated  locally  by  a  quadratic  ray  pencil  with  a  polarization  which 
is  transverse  to  the  incident  ray  direction.  One  notes  that  plane,  cylind¬ 
rical,  conical,  and  spherical  type  incident  wavefronts  are  simply  special 
cases  of  the  arbitrary  quadratic  wavefront.  It  is  further  assumed  in  this 
paper  that  the  field  point  location  and  the  caustics  of  the  incident  ray 
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pencil  are  not  in  the  close  vicinity  of  the  surface,  and  that  the  amplitude 
of  the  incident  ray  optical  field  does  not  exhibit  a  rapid  spatial  variation 
in  the  vicinity  of  the  points  of  reflection  and  diffraction  on  the  convex 
surface.  The  details  of  this  analysis  are  described  in  Sections  II  and  III. 
First,  a  brief  review  of  the  GTD  solution  for  this  general  problem  of  the 
scattering  of  an  arbitrary  ray  optical  electromagnetic  field  by  a  smooth, 
perfectly-conducting  convex  surface  is  presented  in  Section  II  for  the  sake 
of  completeness.  A  uniform  GTD  solution  to  this  problem  is  developed  next 
in  Section  III.  Several  examples  illustrating  the  usefulness  and  accuracy 
of  this  uniform  GTD  solution  are  presented  in  Section  IV. 

II.  THE  GTD  SOLUTION  -  A  REVIEW 

This  review  serves  to  introduce  the  GTD  format  and  notation  that  will 
be  employed  in  the  development  of  the  uniform  GTD  solution. 

The  incident  ray  system  which  strikes  the  convex  surface  produces  a 
system  of  rays  reflected  from  that  surface,  and  the  field  of  these  incident 
and  reflected  rays  constitutes  the  usual  geometrical  optics  (GO)  field. 

The  reflected  ray  merges  with  the  incident  ray  at  grazing  forming  a  shadow 
boundary  which  divides  the  space  exterior  to  the  convex  surface  into  the 
lit  and  shadow  regions,  as  shown  in  Figure  1.  At  grazing,  the  incident  field 
launches  Keller's  surface  rays  which  propagate  along  geodesic  paths  on  the 
convex  surface  while  continually  shedding  energy  via  diffraction  along  the 
forward  tangents  to  the  surface  ray  paths.  The  field  of  these  rays  which 
are  shed  for  diffracted)  from  the  surface  is  known  as  the  surface  diffracted 
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field.  The  GO  field  is  zero  within  the  shadow  region  behind  the  obstacle; 
thus,  the  field  within  the  shadow  zone  is  produced  entirely  by  the  surface 
diffracted  rays.  These  surface  diffracted  rays  may  also  be  present  in  the 
lit  region  if  the  surface  of  the  obstacle  is  closed;  however,  this  field 
is  generally  negligible  compared  to  the  GO  field  if  the  closed  surface  is 
sufficiently  large  in  terms  of  the  wavelength. 

The  GO  solution  for  the  lit  region  is  briefly  reviewed  in  oart  A  of 
this  section.  A  brief  review  of  the  surface  diffracted  ray  solution  for 
the  shadow  zone  is  presented  next  in  part  B.  An  exp(ju)t)  time  dependence 
is  is  assumed  in  this  analysis. 

A.  The  Geometrical  Optics  Field  Solution  for  the  Lit  Region 

Luneberg  [12]  and  Kline  [133  developed  an  asymptotic  high  frequency 
solution  of  Maxwell's  equations  in  which  the  fields  are  expanded  in  inverse 
powers  of  the  angular  frequency,  o).  The  leading  term  in  this  expansion  is 
regarded  as  the  GO  field.  The  details  of  such  an  expansion  are  discussed 
elsewhere  [12,13,14,15]  and  are  only  summarized  here.  According  to  Luneberg 
and  Kline,  the  electric  field  intensity,  T  in  a  source  free,  homogeneous, 
isotropic  medium  can  be  expressed  for  large  u»  by 


n=0  (-jk)' 


fl) 


in  which  r  is  the  position  vector  of  the  field  point;  k  =  w/c;  and  c  is  the 

speed  of  light  in  the  given  medium.  The  coefficients  in  the  above 

expansion  are  determined  by  substituting  Equation  (1)  into  the  vector  Helm¬ 
holtz  equation  satisfied  by  T;  namely,  (V^  +  k^)  T  =  o,  and  by  equating  like 


powers  of  cu.  This  leads  to  the  usual  eikonal  and  transport  equations;  namely, 
|Vi|i|^  =  1  and  +  2v<b*V|  ^-1  “  'respectively. 

The  surfaces  of  constant  iIj  are  referred  to  as  wavefronts,  and  the  family 
of  all  wavefronts  describe  a  system  of  associated  rays  which  are  straight 
lines  in  a  homogeneous  medium.  The  rays  are  everywhere  orthogonal  to  the 
wavefronts  in  an  isotropic  medium.  Integrating  the  transport  equation  for 
the  n=0  case  from  some  reference  point  7^  to  the  field  point  7,  and  expressing 
this  result  in  terms  of  the  Gaussian  curvatures  of  the  wavefronts  at  7^  and 
7  yields  the  GO  field: 


(2) 


The  superscripts  "i"  and  "r"  refer  to  the  incident  and  reflected  GO  fields, 
respectively;  thus,  T’  is  the  incident  GO  field  and  is  the  reflected  GO 
field.  Note  that  and  are  the  principal  radii  of  curvatures  of 

the  incident  or  reflected  wavefront  surface  dA^  at  7^,  and  s^’*^  is  the  distance 
along  the  incident  or  reflected  ray  from  7^  to  7  as  shown  in  Figure  2. 

Next,  requiring  V*f  =  0  leads  to  ($^’’•7^’'"  )  =  0,  which  implies  that  the 

^4  r 

field  in  Equation  (2)  is  polarized  transverse  to  the  ray  directions  (s’)  as 
shown  in  Figure  3.  The  quantity  involving  the  square  root  in  Equation  (2) 

is  the  ray  divergence  factor  which  indicates  the  manner  in  which  the  energy 

spreads  along  the  ray  path;  it  is  a  consequence  of  the  conservation  of  energy 

in  a  ray  tube  (or  pencil).  From  Maxwell's  equation  VxF  =  -  jwyR,  it  follows 

that  the  leading  term  in  the  Luneberg-Kline  expansion  for  the  corresponding 


(A> 


f'^gure  2.  Astigmatic  tube  of  rays 
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magnetic  field  intensity  isTT’’'^~  s’’*"  x  Y^T’**"  in  which  T’’*"  is  as  given 
in  Equation  (’),  and  is  the  characteristic  admittance  of  the  medium. 

From  the  boundary  condition  n  x  =  0  on  the  perfectly-conducting 

ys 

surface,  where  n  is  the  surface  normal  at  the  point  of  reflection  Qp  as  shown 
in  Figure  3,  one  obtains 


T-r 


E'{Qp)  =  E’(Qp)  .  R 


(3) 


Thus,  the  reflected  GO  field  is  given  in  terms  of  the  incident  field  at  Qp 
as 


^■(F)  -r’(Q|,) 


R 


Pi  P2 _ e-jks*^ 

(p1+s’')(p5+s'") 


R  = 


(4) 

f5a;5b) 


where  R  is  the  dyadic  reflection  coefficient,  and  p^,  and  s'"  are  measured 
with  respect  to  the  reference  point  r^  which  is  now  moved  to  Qd.  R  and 
are  the  acoustic  soft  and  hard  reflection  coefficients,  respectively. 

A 

The  unit  vector  e^  in  Equation  (5a)  is  perpendicular  to  the  plane  of  inci¬ 
dence;  whereas,  the  unit  vectors  e’  and  ej|  are  in  the  plane  of  incidence 
as  shown  in  Figure  3.  The  principal  radii  of  curvatures  (PpPg)  of  the  re- 

A  ys 

fleeted  wavefront  and  their  associated  principal  directions  (X^Xp  are  de¬ 
scribed  in  Reference  [16].  The  p!J  and  p^  as  given  in  [16]  are  expressed 
more  compactly  below: 
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- i — .(it 

<■1  ■’(i 
2 


^^—[7,  - 7j  ^~{7,-7j 


gcos2a 


-  sin2«  sin2£i)  cose 


Pg(QR)  ■  "O'  ■  O' 


'ft-y 


4p^(Q Jcos^e'  ■' 
+  1  .  -3-^ - 


1/2 


(6a) 


where 


1  1 
~1  =  7 

Pm 


■  W 


cos^ei 


(6b;6c) 


The  quantities  pj  and  p^  constitute  the  principal  radii  of  curvatures  of  the 

W%  M 

incident  wavefront.  In  Equation  (6a),  P^  and  P^  are  evaluated  at  the  point 
of  reflection  Qj^  on  the  surface.  R-j  and  R^  constitute  the  principal  radii  of 

y\  ^ 

curvatures  of  the  surface  at  Qj^,  and  U|  and  denote  the  corresponding 
principal  surface  directions  at  Qj^.  The  radius  of  curvature  of  the  surface  at 

^  -j  ^  ^ 

Op  is  Pj^;  it  is  measured  in  the  plane  of  incidence  which  contains  s  ,  n  and  t, 

A 

where  t  is  tangent  to  the  surface.  Also,  p^  is  the  radius  of  curvature  of  the 

A  A 

surface  at  Qp  in  the  plane  containing  n  and  the  binormal  vector  b.  The  unit 

A  A  A  A  A 

vectors  t,  n,  b,  and  are  shown  in  Figure  4a  together  with  the  angle 
u)q  between  t  and  U^.  The  unit  vectors  Xj,  X^,  and  the  angle  are  shown 

J  A  A  •  •  A  A 

in  Figure  4b.  The  angle  of  incidence,  0  is  defined  by  n*s  *  -  cos0’  =  -  n*s 
It  is  noted  that  the  GO  representation  of  Equation  (2)  fails  at  caustics 
which  are  the  intersection  of  the  paraxial  rays  (comprising  the  ray  tube 
or  pencil)  at  the  lines  1-2,  and  3-4  as  shown  in  Figure  2.  Upon  crossing 
a  caustic  in  the  direction  of  propagation,  (p  *  +  s  ’jchanqes  sign  under 

the  radical  in  Equation  (2),  and  a  phase  jump  of  +n/2  results  (for  e^^“^ 
time  dependence).  Furthermore,  the  reflected  GO  field,  of  Equation  (4) 
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n  =  UNIT  OUTWARD  NORMAL  VECTOR  TO  THE 
CONVEX  SURFACE  AT  Qr 


Figure  3.  Reflection  by  a  curved  surface. 


/  (a)  (b) 


Figure  4.  Geometry  for  the  description  of  the  wavefront  reflected 

from  a  curved  surface. 


fails  at  and  near  grazing  incidence.  However,  it  is  important  to  note  that 
near  grazing  incidence  (0  -*-ir/2),  and  of  Equation  (6)  approach  the  fol¬ 
lowing  limiting  values. 

(^(Qp)cos9' 


r 

PI 


-o; 


r  1 
P2'^  Pb  » 


(7a; 7b) 


where 


l/p|^  =  (sin^a^/p])  +  (cos^Oq/p^), 


(8) 


.  ^  ys. 

and  p|j  is  the  radius  of  curvature  of  the  incident  wavefront  in  the  (t,b) 
plane  (i.e.,  in  the  plane  tangent  to  the  surface)  at  for  0^-^it/2.  Further- 
more,  the  principal  directions  and  of  the  reflected  wavefront  approach 
the  following  values  for  grazing  incidence: 

X\  (at  0^)  ;  X^  =(-  ^  x  ?f)  ^  b  (at  Q^^)  .  {9a;9b) 

The  limiting  values  in  Equation  (9a;b)  are  independent  of  o^.  The  total 
GO  field,  E  at  Pj^  (see  Figure  1)  in  the  lit  region  is  the  sum  of  the  incident 
and  reflected  GO  ray  fields;  hence. 


E(Pl)  -^'(Pl)  +  IVQr) 


r  r 
Pi  P2 


<1  (p!;+s»')(p5+s^) 


(10) 


B.  The  Surface  Diffracted  Ray  Field 
Solution  for  the  Shadow  Region 

The  incident  ray  at  grazing  launches  a  set  of  surface  rays  which  propa¬ 


gate  along  a  geodesic  path  on  the  convex  surface  thereby  carrying  energy 
into  the  shadow  region.  The  field  associated  with  these  surface  rays  attenuates 
(i.e.,  decays  exponential Iv)  as  they  propagate  due  to  a  continuous  shedding 


m 


(diffraction)  of  rays  along  the  forward  tangents  to  the  geodesic  ray  path, 
as  shown  in  Figure  5.  An  analysis  of  this  surface  diffracted  field  is  dis¬ 
cussed  in  detail  elsewhere  [1,2,15];  hence,  only  the  essential  features  are 

A  A 

summarized  here.  Let  n^  and  t-j  denote  the  unit  normal  and  tangent  vectors 
at  the  point  of  grazing  incidence,  Q^;  likewise,  let  n^  and  t^  denote  the 
unit  normal  and  tangent  vectors,  respectively,  at  the  point  of  tangential 
shedding  of  the  surface  diffracted  ray  at  Q^.  The  diffracted  ray  pencil 
possesses  caustics  at  and  P^.  as  seen  in  Figure  5;  these  caustics  are  the 
same  as  those  at  3-4  and  1-2  of  the  typical  ray  pencil  in  Figure  2.  The 
diffracted  electric  field  arriving  at  may  be  represented  as  a  ray  op¬ 
tical  field;  hence  is  given  by  Equation  (2)  with  the  superscripts  "i,r" 
replaced  by  "d"  to  denote  the  "diffracted"  ray.  Then  and  become  the 
caustic  distances  associated  with  the  diffracted  wavefront  at  some  reference 
point  (r^)  which  lies  between  and  P^,  as  in  Figure  2.  However,  if  the 
reference  position  r^  is  moved  back  to  the  caustic  at  Q^,  then  p^  -►  0,  and 
now  becomes  the  distance  from  to  P^.  Let 

j  _  -jk\l)'^{r  )  r—r  . 

lim  r(r^)e  °  Jp?  =  l’(Q^).T{Qi,Q2)  .  (11) 


Therefore,  the  electric  field  fr^)  at  a  point  P^  in  the  shadow  region  (see 
Figure  1)  becomes 


^(Pj)  ~  ■P’(Qt)  •  t(QTQ2) 


d 


(12a) 


is  a  "dyadic  transfer  function"  which  relates  the  field  diffracted 
from  to  the  field  incident  at  Q^.  This  dyadic  quantity  is  expressed  as[3,15j 


12 


Figure  5.  Diffraction  by  a  smooth  curved  surface 
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TfQijQ^)  “  ^1^2  ^  "^h  * 


(12b) 


in  which 


Q2  5 

-jkt-  /  ap(t')df  - — 

fi  I  Ql  dn(Qi) 

h  L 


dn(Qi)  I 
dn(Q2) 


and  t  is  the  geodesic  arc  length  from  to  Q2  on  the  surface.  It  is  clear 
from  Equation  (12b)  that  the  ray  field  in  Equation  {12a)  is  polarized  transverse 
to  the  ray  path.  It  is  assumed  in  Equations  (12a;b)  that  the  02  and  b^  com¬ 
ponents  of  propagate  independently  of  each  other.  The  sum  in  Equation 
(12c)  indicates  that  the  surface  ray  field  is  actually  composed  of  a  set 
of  surface  ray  modes  as  indicated  earlier,  and  p  refers  to  the  modal  index. 

The  superscripts  s  and  h  in  Equations  (12b;c)  d|note  the  acoustic  soft  and 

u 

hard  type  field  contributions,  respectively.  OpfO-j)  are  referred  to  as  the  ' 

surface  dif-fr^ction  coefficients  which  describe  t^e  diffraction  at  Q-j. 

h  h  ^ 

The  forms  of  DpfO^)  and  Dp(Q2)  must  be  identical  via  reciprocity.  The  factor 


^2  ! 


dn(Qi) 


-jkt  -  /a5(t’)df 


is  the  ratio  of  the  surface  ray  field  incident  at  Q2  (prior  to  diffraction 
from  Qp)  to  the  surface  ray  field  launched  by  the  incident  field  at  Q^. 
This  factor  is  obtained  by  integrating 


^  (ag  dn)  = 
h 


s 

2(ap  dn)  , 
h 
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from  to  Q^,  with  being  the  intensity  of  the  p  soft  and  hard  surface 
ray  field  and  dt  is  tfle  incremental  arc  length  between  and  Q2.  This  ex¬ 
pression  is  based  on  the  assumption  that  energy  is  conserved  in  the  surface 
ray  pencil  between  and  Q^,  and^that  the  surface  ray  decays  exponentially 
with  an  attenuation  coefficient  0^.  The  factor  denotes  the  dominant 

phase  delay  of  the  surface  ray  field  from  Qi  to  Q^. 

The  form  of  the  solution  in  Equations  (12a;b;c)  has  been  verified  via 

asymptotic  solutions  to  appropriate  canonical  problems  [3,17];  a^so,  ^ 

h  h 

these  canonical  solutions  lead  to  the  specific  expressions  for  Dp  and  0^ 
which  to  first  order  are  given  by  [3,17] 


-j^2 

e _ 


h 


(13a;13b) 

(14a;14b) 


The  Ai  and  Ai'  denote  the  Keller  type  Airy  function  and  it's  derivative 
[3,15]  respectively.  The  values  q  and  q  are  those  for  which  Ai(-q  )  =  0 

K  K  K 

and  Ai'(-qp)  =  0.  The  values  of  Al(-qp)  and  Ai'(-qp)  are  given  in  [9,18], 

It  is  noted  that  refers  to  the  surface  radius  of  curvature  in  the  ray 

A  A 

direction  (i.e.,  in  the  n,t  plane).  In  Equation  (12c),  the  inclusion  of 
only  the  first  couple  of  modes  (i.e.,  p=l,2)  is  sufficient  for  obtaining 
accurate  results  in  the  deep  shadow  region.  However,  the  result  in  Equations 
n2a;h;cl  fails  at  and  near  the  shadow  boundary;  it  also  fails  near  the 
caustic  at  Qp. 
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III.  THE  CONSTRUCTION  OF  A  UNIFORM  6T0  SOLUTION 


The  GTD  solution  of  Section  II  fails  at  and  near  the  shadow  boundary 
because  the  GTD  ray  optical  field  description  is  not  valid  there.  Conse¬ 
quently,  one  must  employ  a  uniform  asymptotic  field  approximation  which  remains 
valid  within  the  shadow  boundary  transition  region.  A  uniform  theory  must 
therefore  basically  depart  from  the  pure  ray  optical  field  approximation 
inherent  in  the  GTD  in  order  to  correct  for  the  failure  of  the  GTD  within 
the  shadow  boundary  (SB)  transition  region;  whereas,  away  from  this  SB  tran¬ 
sition  region,  it  must  reduce  to  the  usual  GTD  solution  where  the  latter 
is  indeed  valid.  The  precise  manner  in  which  a  uniform  theory  accomplishes 
such  a  task  differs  with  each  ansatz.  While  different  formulations  of  uni¬ 
form  solutions  might  even  lead  to  the  same  answers  exactly  at  the  SB,  and 
also  exterior  to  the  SB  transition  region  where  they  all  must  reduce  to  the 
GTD  solution,  their  behavior  within  the  transition  region  may  not  necessarily 
be  the  same. 

As  mentioned  in  Section  I,  the  ansatz  employed  in  this  paper  is  based 
on  an  asymptotic  solution  given  by  Pathak  [4]  for  the  canonical  problem  of 
plane  wave  scattering  by  a  smooth,  perfectly-conducting  convex  cylinder. 

This  ansatz  leads  to  a  uniform  GTD  solution  for  the  general  problem  of  the 
scattering  of  a  ray  optical  electromagnetic  field  by  a  smooth,  perfectly- 
conductina  convex  surface  of  any  shape,  such  that  the  solution  thus  obtained 
is  convenient  and  accurate  for  engineering  applications.  The  starting  point 
in  the  development  of  this  general  solution  is  the  uniform  result  given  in 
Equations  (A-lfia;b)  of  the  Appendix,  for  the  far  zone  fields  of  a  scalar 
point  source  radiating  in  the  presence  of  an  acoustic  soft  or  hard,  smooth 
convex  cylinder.  The  result  in  Equations  {A-16a;b)  is  developed  in  the  Ap¬ 
pendix  from  the  uniform  asymptotic  solution  of  the  canonical  problem  of 
[4], the  latter  solution  is  also  summarized  in  Equations  (A-la;b)  of  the  Ap¬ 
pendix.  The  reader  is  referred  to  the  Appendix  for  details. 
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First,  it  is  observed  that  the  uniform  result  in  Equations  (A-16a;b) 
is  already  in  the  convenient  ray  type  format;  however,  the  field  itself  is 
not  a  ray  optical  field  within  the  SB  transition  region.  The  scattered  field 
in  Equations  (A-lfia;b)  propagates  along  the  GO  reflected  ray  path  in  the 
lit  region  (which  also  includes  the  lit  portion  of  the  SB  transition  region!, 
and  it  propagates  along  the  surface  diffracted  ray  path  in  the  shadow  region 
(including  the  shadowed  portion  of  the  SB  transition  region).  This  uniform 
result  properly  reduces  to  the  GO  soft  or  hard  type  ray  field  solution  in 
the  lit  region  which  lies  outside  the  SB  transition  region,  and  also  to  the 
soft  or  hard  type  surface  diffracted  ray  field  solution  in  the  shadow  region 
which  lies  exterior  to  the  SB  transition  region.  The  preceding  reduction 

A 

to  the  GTD  solution  is  easily  verified  by  noting  that  the  F  and  functions 

u 

which  occur  in  Equations  (A-16a;b)  of  the  Appendix  take  the  following  lim¬ 
iting  values,  when  the  field  point  moves  exterior  to  the  SB  transition  region: 


(which  is  true  both,  in  the  lit  and 
shadow  zones  exterior  to  the  SB 
transition  region. 

(15) 


(which  is  true  in  the  lit  zone\ 
exterior  to  the  SB  transition  1 
region.  / 


(16a) 


N 


P3(«) 


ilL  Xr,  a 

e  e  ^ 


f  P=1  2[Ai'(-qp)]' 


N  fi  ” 

e _  r  e  e  ^  _ 

P=1  2  qp[Aif-qp)]^ 


Phf6) 


which  is  true  in 
for  5»0  the  shadow  zone 

exterior  to  the  SB 
\ transition  region. 


(16b) 


The  N  in  the  summation  of  Equation  (16b)  is  identical  to  the  N  in  Equation 
(^2c)\  it  is  noted  that  N=2  generally  provides  sufficient  accuracy  when  6»0. 
Furthermore,  the  limit  of  the  result  in  Equation  fA-16a)  as  the  field  point 
approaches  SB  from  the  lit  side  is  identical  to  the  limit  of  the  result  in 
Equation  (A-16b)  as  the  field  point  approaches  SB  from  the  shadow  side;  thus, 
the  total  field  is  continuous  across  SB.  It  is  noted  that  as  the  field  point 
approaches  SB: 

j|l  j(|  +  a) 

F(o)  IjiTal  -  2o  e  ^je  ^  ,  for  a-K)  (ct=0  on  SB).  (17) 

Also,  the  limiting  value  of  the  field  at  the  SB  is  more  easily  evaluated 

A 

if  one  defines  2^(6)  in  terms  of  the  related  functions  p*(6)  and  q*(6)  as 
in  [4]: 


PJ6)  = 


P*f6) 

q*(«) 


-il  “J4 
.  ^  e  ^ 


(Note  that  6=0  at  SB).  (18) 


From  the  above  limiting  forms  of  F(o)  and  2^(6),  one  notes  that  the  F(o) 

h 

term  in  Equations  (A-16a;b)  plays  a  dominant  role  in  the  immediate  neighbor¬ 


hood  of  SB,  and  it  is  entirely  responsible  for  ensuring  the  continuity  of 
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the  total  field  at  SB.  On  the  other  hand  the  P5f6)  term  in  Equations  (A- 

u 

16a;b)  plays  a  dominant  role  as  the  field  point  moves  far  from  SB  (since 
|l-Ffol|  -*■  0),  and  it  is  therefore  entirely  responsible  for  reducing  Equations 
fA-lfia;b)  uniformly  to  the  GTD  solution,  exterior  to  the  SB  transition  region, 
for  the  case  of  the  far  zone  field  of  a  scalar  point  source  radiating  in 
the  presence  of  an  acoustic  soft  or  hard  convex  cylinder.  A  more  complete 
discussion  on  the  role  of  the  functions  F  and  is  given  in  Reference  [4], 
in  connection  with  the  development  of  a  uniform"asymptotic  solution  to  the 
canonical  problem  of  plane  wave  scattering  by  a  convex  cylinder.  It  is  im¬ 
portant  to  note  that  the  F  function  involves  a  Fresnel  integral  which  is 
well  tabulated;  and  P^  is  a  Fock  type  integral  (involving  Airy  functions) 
which  is  also  tabul  at&d. 

It  is  observed  that  the  GTD  solution  of  Section  II,  which  is  a  first 
order  asymptotic  solution  to  terms  in  inverse  powers  of  kp^,  is  valid  for 
cylindrical,  spherical,  or  any  other  smooth  convex  shape.  It  is  also  valid 
for  torsional  surface  rays  in  that  effects  of  torsion  do  not  occur  explicitly 
to  first  order  for  the  scattering  problem  considered  here.  In  addition, 
it  is  observed  that  a  simple  relationship  exists  in  the  GTD  solution  between 
the  vector  electromagnetic  and  the  scalar  acoustic  problems  to  the  given 
order  of  approximation;  namely,  R  and  f  of  the  GTD  solution  in  Equations 
{5a;b)  and  (12a;b)  are  a  simple  combination  of  the  corresponding  scalar  or 
acoustic  (soft  and  hard)  functions  R^,  and  T^  for  the  lit  and  shadow  regions, 
respectively.  The  last  two  observations  in  l^egard  to  the  GTD  solution  in¬ 
cluding  it's  validity  for  torsional  surface  rays  have  been  verified  via  the 
rigorous  asymptotic  solutions  to  several  canonical  problems  [3,18]  .  Based 
on  these  observations  and  the  previous  observation  that  the  uniform  asymp¬ 
totic  result  in  Equations  (A-16a;b)  of  the  Appendix  dealing  with  the  far 


zone  field  of  a  point  source  radiating  in  the  presence  of  an  acoustic  soft 
or  hard  convex  cylinder  is  already  in  a  ray  format,  it  is  reasonable  to  con¬ 
jecture  the  following  in  regard  to  the  uniform  GTD  solution: 

(a)  The  uniform  GTD  solution  for  the  general  electromagnetic  problem 
of  the  scattering  of  a  ray  optical  electromagnetic  field  incident  on 
a  smooth,  perfectly-conducting  arbitrary  convex  surface  can  also  be 
expressed  in  a  ray  type  format  like  the  ordinary  GTD  solution  of  Sec¬ 
tion  II  (see  Equations  (10)  and  (12a)). 

(b)  The  dyadics  R  and  f  in  Equations  (5a;b)  and  (12a;b)  of  the  or¬ 
dinary  GTD  solution  may  be  replaced  by  the  more  general  I  and  f  dyadics 
for  the  uniform  GTD  solution.  Of  course,  R  -*■  R  must  be  true  in  the 
lit  region  outside  the  SB  transition  region  and  R  f  must  be  true 

in  the  shadow  region  outside  the  SB  transition  region.  Furthermore, 
the  functional  forms  of  R  and  f  are  assumed  to  be  the  same  for  cyl¬ 
indrical,  spherical,  or  anv  other  convex  shape,  as  is  true  of  the  R 
and  t  dyads  in  the  GTD  solution. 

(c)  The  uniform  GTD  solution  for  the  electromagnetic  case  may  also 

be  simply  expressed  in  terms  of  the  corresponding  scalar  (or  acoustic) 
soft  and  hard  cases,  respectively,  as  in  the  GTD  solution  [  see  Equa¬ 
tions  (5a) ,  and  (12b)]  . 

Thus,  the  form  of  the  uniform  GTD  solution  for  the  total  electric 
field,  E  may  be  expressed  via  the  con.i'ctures  in  (a)  and  (b)  above,  as: 
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Eor  P^  in  the 
shadow  region.  (19b) 
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The  incident  ray  field  is  the  same  as  in  Equation  (2)  or  (10).  From  con¬ 
jecture  (c)  above,  one  may  express  R  and  f  in  terms  of  their  corresponding 
acoustic  soft  and  hard  functions  as: 


R  =  Rg  e^  +  R^ 


(20) 


T 


+  Tf,  n^n^ 


(21) 


The  results  in  Equations  (19a;b)  together  with  the  dyadics  in  Equations  (20) 

and  (21)  are  expressed  in  the  same  incident,  reflected  and  surface  diffracted 

ray  coordinates  as  those  employed  in  the  GTD  solution  of  Section  II.  The 

subscripts  s  and  h  in  Equations  (20)  and  (21)  refer  to  the  acoustic  soft 

and  hard  cases,  as  before.  From  Equations  (20)  and  (21),  and  conjecture 

(b)  above,  an  obvious  choice  of  R  and  T  is  directly  available  from  the 

h  h 

acoustic  soft  or  hard  solution  of  Equations  (A-16a;b)  in  the  Appendix.  Thus, 
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for  the  lit 
region 
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■jkt 


for  the 
shadow  region. 

(23) 

The  various  parameters  occurring  in  Equations  (22)  and  (23)  are  thus  defined 
below  as  in  Equations  (A-lTa)  through  (A-llj.)  of  the  Appendix. 
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For  an  electromagnetic  spherical  wave  (or  point  source)  type  illumination, 

r  .  r  _ ,  _d  d  . 


(19a;b)  such  that 


and  for  the  field  point  in  the  far  zone  where  s  »p-,  ^  and  s  »po  in  Equations 
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,  the 


"distance"  parameters  L ^  and t °  in  Equations  (24f;24g)  are  given  for  this 
special  case  by  f  ^  =  s'"  and  1- =  s^  as  in  Equations  (A-16a;b).  It  is  noted 
that  s'",  s*^  andp^  appearing  in  Equations  (19a;b)  correspond  to  j, ,  s^  and 
in  Equations  (A-16a;b),  respectively.  Therefore,  in  order  to  complete 
the  solution  in  Equations  (19a;b)  for  the  general  case,  one  need  only  specify 
the  parameters  and  To  recapitulate,  the  general  case  deals  with  an 
arbitrary  incident  ray  optical  electromagnetic  field,  an  arbitrary  convex 


surface,  and  a  near  zone  field  point  (for  which  the  field  point  may*^ be  sev¬ 
eral  wavelengths  from  the  scattering  surface,  but  yet  not  far  in  comparison 
with  the  characteristic  dimensions  of  the  scatterer).  As  mentioned  previ¬ 
ously  in  Section  I,  the  plane,  cylindrical,  conical  and  spherical  wavefronts 
are  special  cases  of  the  arbitrary  quadratic  wavefront  approximation  implied 
in  the  ray  optical  description  of  the  incident  field.  Recalling  that 
the  role  of  the  F(X^)  and  F(X^)  functions  is  to  ensure  the  continuity 
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of  the  total  field  at  SB,  one  may  then  evaluate  and  X*^  appearing  in 
X*"  and  X^,  respectively,  by  actually  enforcing  the  continuity  of  the 
total  field  at  SB.  This  procedure  is  exactly  analogous  to  that  employed 
previously  by  Kouyoumjian  and  Pathak  jlS]  in  their  development  of  a  uniform 
GTD  solution  for  the  diffraction  of  an  arbitrary  ray  optical  electromagnetic 
field  by  an  edge  in  an  otherwise  smooth  surface.  Let  P^g  denote  a  field 
point  on  SB.  The  continuity  of  the  total  field  at  SB  requires  that 
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in  which  E^Pj^)  and  E(P^)  are  given  by  Equations  (19a)  and  (19b}.  Employing 


the  limiting  form  of  Equation  (17)  for  the  F  function,  and  the  definition 

le  Pg  fun( 

allows  one  to  write  Equation  (25)  as: 


given  in  Equation  (18)  for  the  P^  function,  into  Equations  (22)  and  (23), 
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With  the  distance  from  Q-j  to  P^g  defined  as  s,  one  obtains. 
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It  is  noted  that  if  PL^P5g,  then  P^-^O  (see  Equatii 

(7a)),  e„>n^,  e„+np  and  e^+b^.  Likewise,  if  P5-^P5g,  then  5  +0,  Q^-^-Qp 
t>-0,  n^^np  and  b^+b^.  Furthermore,  at  SB,  02=02*  l^be  functions  p*(6) 
and  q*((S)  are  continuous  everywhere  including  (S=0.  Let  the  incident  field 
be  described  by  a  diverging  wavefront  in  the  direction  of  propagation  at 
(i.e.,  at  grazing);  thus,  one  may  write, 

-jks„ 
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Pi  P2 


(Pl%)(P2+s^) 
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where  s^  denotes  the  distance  s^  from  some  reference  point  at  r^  to  the  point 

_  A>.  Ak 

O-j.  The  quantities  3*(r^)  and  3„(r^)  which  are  the  b-j  and  n-j  directed  ampli¬ 


tudes  of  E  at  the  reference  location  (r^)  are  assumed  known.  It  follows 


from  Equation  (28)  that 
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The  incident  wave  caustic  distances  o’  ^  are  measured  from  the  reference 


point  at  r^  to  the  respective  caustic  locations.  Incorporating  Equations 
(28)  and  (29)  into  (26)  yields 


(pj+So)(pJ+So) 


s(p2+s) 


SB  'SB  (p;+[So+s])(pJ+[s^+s])  pi 


at  SB. 


(30) 


Since  the  distance  parameter  is  a  slowly  varying  quantity  near  SB,  and 

since  (1-F)  vanishes  sufficiently  rapidly  as  the  field  point  moves 

far  from  SB;  it  is  convenient  to  use  the  value  in  Equation  (30)  for  L*-  and 


i  even  away  from  SB.  It  is  assumed  of  course  in  the  present  development 


(as  pointed  out  in  Section  I),  that  the  field  point  location  and  the  caustics 
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of  the  incident  wavefront  are  not  in  the  close  vicinity  of  the  surface,  and 
that  the  field  point  itself  is  not  in  the  neighborhood  of  any  caustics  as¬ 
sociated  with  the  incident  and  scattered  rays.  Furthermore,  it  is  assumed 
that  the  amplitude  of  the  incident  ray  optical  field  does  not  exhibit  a  rapid 
spatial  variation  in  the  vicinity  of  the  points  of  reflection  and  diffraction 
on  the  surface.  This  completes  the  construction  of  the  uniform  GTD  solution 
in  Equations  (19a;b)  for  a  diverging  wavefront. 

If  the  incident  wavefront  is  of  the  converging,  or  converging-diverging 
type,  then  the  parameters  in  Equation  (30)  can  become  negative.  It 
has  not  been  investigated  in  detail  how  the  general  solution  can  be  completed 
when  becomes  negative.  However,  if  one  of  the  principal  directions 
of  the  incident  wavefront  coincides  with  one  of  the  principal  planes  of  the 
surface  at  grazing,  then  one  can  treat  a  converging  or  converging-diverqinq 
type  wavefront  for  which  by  replacing  F(X*'’^)  with  F*(lx^’^l).  The 

*  on  F*f  jx**’^!)  denotes  the  complex  conjugate  operator.  The  use  of  F*(  Ix^’^l) 
when  L*"’^<0  leads  to  a  continuous  total  field  at  SB  in  this  case,  and  it's 
use  may  be  justified  as  in  the  edge  diffraction  problem  [19]  via  an  analytic 
continuation  procedure  to  include  negative  values  of  (or  while 

simultaneously  satisfying  the  radiation  condition  for  the  scattered  field. 

Finally,  exterior  to  the  SB  transition  region,  the  uniform  GTD  result 
of  Equations  (19a;b)  does  indeed  recover  the  ordinary  GTD  result  of  Equations 
(10)  and  (12a)  as  may  be  verified  by  employing  the  limiting  forms  of  Equations 
(15)  and  (16a;b)  into  Equations  (22)  and  (23). 

i 

i 

i 
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It  is  noted  that  a  solution  for  the  SB  transition  region  pertaining 
to  the  same  general  problem  as  the  one  considered  here  has  also  been  given 
by  G.  James*;  however,  his  solution  which  is  presented  without  sufficient 
details  is  not  uniform  in  that  it  does  not  properly  reduce  to  the  GTD  solu¬ 
tion  exterior  to  the  SB  transition  region.  Furthermore,  that  solution  em¬ 
ploys  a  "pseudo"  ray  path  for  the  reflected  field  in  the  lit  region;  such 
a  path  does  not  satisfy  the  generalized  Fermat’s  principle.  On  the  other 
hand,  as  is  well  known,  the  GO  reflected  ray  path  which  is  employed  in  this 
paper  indeed  satisfies  the  generalized  Fermat's  principle. 

IV.  DISCUSSION  AND  NUMERICAL  RESULTS 

A  uniform  GTD  solution  has  been  obtained  for  the  problem  of  the  scat¬ 
tering  of  a  ray  optical  electromagnetic  field  by  a  smooth,  perfectly-conducting 
convex  surface  as  shown  in  Figure  1.  This  result  is  explicitly  given  in 
Equation  {19a)  for  the  lit  region,  and  in  Equation  (19b)  for  the  shadow  region, 
together  with  the  parameters  in  Equations  (?0)  through  {24g),  and  also  Eq¬ 
uation  ITOI.  While  the  behavior  across  only  a  single  shadow  boundary  (SB) 
is  discussed  in  this  paper  in  connection  with  the  open  convex  surface  of 
Figure  1,  the  present  theory  can  just  as  easily  treat  the  scattering  by  a 
closed  convex  surface.  Basically,  one  treats  a  closed  convex  surface  via 
the  uniform  GTD  in  the  same  manner  as  one  would  via  the  ordinary  GTD;  thus, 
the  only  difference  between  the  two  approaches  is  that  the  R  and  f  dyads 
in  the  GTD  solution  are  replaced  by  the  more  general  and  more  accurate  ^ 
and  t  dyads  in  the  uniform  GTD  solution.  It  is  noted  of  course  that  the 


*G.L.  James,  Geometrical  Theory  of  Diffraction  for  Electromagnetic  Wayes, 


published  by  Peter  Peregrinus  Ltd.,  Southgate  House,  Steyenage,  Herts.  SGI 
IHQ,  England,  1976. 
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usual  GTD  solution  fails  within  the  shadow  boundary  transition  regions;  where¬ 
as,  the  uniform  GTD  solution  does  not.  When  a  surface  of  revolution  is  il¬ 
luminated  by  a  plane  wave  which  is  incident  along  the  axis  of  revolution 
of  the  surface,  a  caustic  of  the  surface  diffracted  rays  lies  on  this  axis. 

The  uniform  GTD  solution  fails  (as  does  the  usual  GTD  solution)  in  the  neigh¬ 
borhood  of  such  a  caustic.  However,  if  the  field  point  is  in  the  near  zone 
such  that  the  caustic  and  the  shadow  boundary  directions  are  widely  sepa¬ 
rated,  then  one  can  employ  the  method  of  equivalent  ring  currents  to  evaluate 
the  field  in  the  neighborhood  of  such  a  caustic.  The  equivalent  currents 
in  this  case  are  found  indirectly  from  the  uniform  GTD  solution.  Such  a 
procedure  will  be  reported  in  a  separate  paper  together  with  another  approxi¬ 
mate  technique  which  would  yield  the  field  near  caustics  for  the  special 
case  when  the  caustic  and  the  shadow  boundary  directions  tend  to  coincide. 

The  latter  special  case  arises  when  the  field  point  is  near  the  axis  but 
in  the  far  zone  behind  the  surface  of  revolution.  The  other  restrictions 
on  the  uniform  GTD  solution  are  mentioned  at  the  end  of  both.  Sections  I 
and  III,  respectively. 

It  is  interesting  to  note  that  the  ^  of  Equation  (24a)  may  be  approxi¬ 
mated  for  convenience  by  -2  m(Qp)cos8’  (upon  arbitrarily  setting  f"^^^(Qj^)=l) 
without  affecting  the  accuracy  of  the  solution.  Furthermore,  it  is 
easily  verified  that  the  parameter  of  Equation  (30)  simplifies  to  give 
the  following: 
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s'+s 


for  a  2-0  cylindrical  wave  M.e.,  line  source!, 
or  an  obliquely  incident  3-D  spherical  wave 
(i.e.,  point  source)  illumination. 


s  ,  for  a  plane  wave  illumination  (s^®).  (31) 

In  the  above  expression,  s'  denotes  the  distance  from  the  source  point  to 
the  point  of  grazing  incidence  at  Q-j.  For  the  cylindrical  wave  illumination 
s'  and  s  are  to  be  interpreted  as  the  distances  in  the  complete  2-0  problem. 
With  L*"’^  as  in  Equation  (31)  for  plane,  cylindrical,  or  spherical  wave  il¬ 
lumination,  it  is  easily  verified  that  the  uniform  GTD  result  of  Equations 
(19a;b)  satisfies  reciprocity.  Finally,  this  uniform  GTD  solution  is  simple 

and  accurate  to  use  since  it  is  given  in  terms  of  the  F  and  type  functions 

h 

which  are  tabulated.  These  important  aspects  of  the  present  uniform  GTD 
solution  are  illustrated  below  by  applying  it  to  several  interesting  prob¬ 
lems. 

Consider  an  infinitely  long  circular  cylinder  of  radius  (a)  as  shown 
in  Figure  6a.  The  geometrical  optics  field  for  this  geometry  is  shown  in 
Figure  6b  for  an  electric  dipole  (acoustic  soft  boundary  condition)  and  in 
Figure  6c  for  a  magnetic  dipole  (acoustic  hard  boundary  condition)  mounted 
parallel  to  the  cylinder  axis.  As  is  well  known,  the  geometrical  optics 
field  breaks  down  near  the  shadow  boundaries.  This  is  readily  apparent  in 
these  two  figures.  Nevertheless,  this  solution  has  often  been  applied  to 
obtain  a  result  for  the  scattering  from  a  cylinder.  The  uniform  GTD  solu¬ 
tion  presented  in  this  paper,  however,  can  be  quickly  and  accurately  computed. 
The  uniform  GTD  result  is  compared  with  the  geometrical  optics  result  in 
Figures  6b  and  c.  Obviously,  a  much  more  complete  result  for  the  problem 


Figure  fib.  Comparison  of  the  geometrical  optics  and  uniform  geometrical 
theory  of  diffraction  radiation  patterns  for  an  electric  dipole 
in  the  presence  of  an  infinitely  long  circular  cylinder  with 

a  =  lx  and  p'  =  3x. 
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Figure  6c.  Comparison  of  the  geometrical  optics  and  uniform  geometrical 
theory  of  diffraction  radiation  patterns  for  a  magnetic  dipole  in 
the  presence  of  an  infinitely  long  circular  cylinder 
with  a  =  1  ^  and  p  '  =  2X  . 
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in  Figure  6a  is  obtained  using  the  uniform  GTD  presented  here.  The  validity 
of  this  result  is  shown  by  comparing  it  with  the  exact  eigenfunction  solution. 
The  exact  and  the  uniform  GTD  results  for  the  two  polarizations  are  compared 
in  Figure  7.  The  agreement  is  excellent  and  hence  confirms  the  validity 
of  the  uniform  GTD  solution. 

The  next  example  is  of  an  electric  dipole  in  the  presence  of  a  finite 
circular  cylinder,  as  illustrated  in  Figure  8a.  Measurements  have  been  made 
on  this  satellite  shape  by  Bach  [20]  and  are  used  for  comparisons  with  the 
calculated  results  in  Figure  8b.  The  pattern  is  taken  in  the  x-y  plane. 

The  end  caps  of  the  cylinder  do  not  have  a  large  effect  in  this  plane  and, 
therefore,  do  not  need  to  be  included.  Note  that  the  dipole  is  aligned 
parallel  to  the  x  axis.  This  causes  the  cylinder  to  be  illuminated  by  a 
slowly  varying  field.  The  measured  and  calculated  results,  however,  are 
in  good  agreement  within  experimental  accuracy. 

An  example  of  a  magnetic  dipole  in  the  presence  of  an  elliptic  cylinder 
is  considered  next,  as  illustrated  in  Figure  9a.  The  pattern  is  a  conic 
cut  about  the  cylinder  axis.  The  dipole  is  parallel  to  the  cylinder  axis. 

This  could  represent  a  slot  mounted  near  an  aircraft  fuselage,  engine,  or 
store.  The  calculated  result  is  compared  against  a  result  obtained  from 
a  moment  method  solution*  in  Figure  9b.  The  agreement  is  very  good,  again, 
verifying  the  validity  of  the  present  solution. 

*This  solution  has  been  kindly  furnished  by  Dr.  N.  Wang. 


Figure  7a.  Comparison  of  the  exact  eigenfunction  and  uniform  geometrical 
theory  of  diffraction  radiation  patterns  for  an  electric  dipole  in 
the  presence  of  an  infinitely  long  circular  cylinder 
with  a  =  and  p  '  =  2x. 
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Figure  7b.  Comparison  of  the  exact  eigenfunction  and  uniform  geometrical 
theory  of  diffraction  radiation  pattern  for  a  magnetic  dipole  in 
the  presence  of  an  infinitely  long  circular  cylinder 
with  a  =  lx  and  p  '  =  2x . 
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Figure  qb  Comparison  of  a  moment  method  and  uniform  geometrical  theory 
of  diffraction  radiation  pattern  for  a  magnetic  dipole  in  the 
presence  of  an  infinitely  long  circular  cylinder.  The  pattern 
IS  a  conic  cut  about  the  cylinder  aixs. 
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In  order  to  validate  this  solution  in  terms  of  a  more  complex  surface 
consider  the  circular  cone  configuration  illustrated  in  Figure  10.  A  half¬ 
wave  dipole  is  located  in  the  near  zone  of  the  cone,  and  the  near  zone  field 
is  measured  as  the  receiver  moves  azimuthly  around  the  axis  of  the  cone. 

Both  a  vertical  and  horizontal  dipole  are  treated  with  the  resulting  patterns 
shown  in  Figures  11a  and  11b,  respectively.  In  each  case  the  calculated 
and  measured  results  are  in  good  agreement.  Note  that  the  receiver  polari¬ 
zation  was  aligned  with  that  of  the  transmitter  in  both  cases. 

The  last  example  is  used  to  show  a  cylinder  being  illuminated  by  a 
complex  wavefront  that  can  be  represented  in  terms  of  an  astigmatic  tube 
of  rays.  The  source  of  the  astigmatic  tube  of  rays  or  the  quadratic  ray 
pencil  which  impinges  on  the  cylinder  is  the  edge  diffracted  field  of  a  plate 
mounted  on  the  cylinder  such  that  the  cylinder  is  not  in  the  shadow  boundary 
transition  region  of  the  edge  diffracted  fields.  The  geometry  used  is  il¬ 
lustrated  in  Figure  12.  The  antenna  is  a  slot  mounted  parallel  to  the  cylinder 
axis  in  the  center  of  the  plate.  The  calculated  and  measured  results  for 
this  configuration  are  shown  in  Figure  12.  There  is  good  agreement  between 
measured  and  calculated  results,  thus  confirming  the  validity  of  the  uniform 
GTD  solution  for  a  cylinder  illuminated  by  a  more  general  wavefront. 
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Figure  11a,  Comparison  of  measured  fdashed  curve)  and  calculated  (solid 
curve)  radiation  patterns  for  an  electric  dipole  mounted  parallel 
to  the  z  axis  of  a  circular  cone. 
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Figure  11b.  Comparison  of  measured  (dashed  curve)  and  calculated  (solid 
curve)  radiation  patterns  for  an  electric  dipole  mounted  paralle 
to  the  y  axis  of  a  circular  cone. 


Figure  1?.  Comparison  of  measured  and  calculated  radiation  patterns  for 
a  slot  mounted  in  a  plate-cylinder  configuration  with  the 
slot  parallel  to  the  cylinder  axis. 
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APPENDIX 


FAR  ZONE  FIELDS  OF  A  SCALAR  POINT  SOURCE  RADIATING  IN 
THE  PRESENCE  OF  A  SMOOTH  CONVEX  CYLINDER 


An  approximate  uniform  asymptotic  result  for  the  field  (u^)  exterior 
to  a  smooth  convex  cylinder  illuminated  by  an  incident  plane  wave  (u^)  is 
presented  in  Equations  (45a)  and  (44a)  of  [4].  It  is  repeated  here  as  fol¬ 
lows: 


U2(P)~U2(P)+U2(Qr) 
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for  P  in  the  shadow  region.  (A-lb) 
The  subscripts  "s"  and  "h"  refer  to  the  acoustic  soft  and  hard  cases, 
respectively;  whereas,  the  subscript  2  refers  to  the  two-dimensional  (2-D) 
nature  of  the  problem.  The  points  P,  Qj^,  and  Qj^  in  Equations  (A-la;b) 


are  illustrated  in  Figure  A-1  for  the  case  when  0^  =  tt/2;  i.e.,  for  "normal 


incidence".  Also,  the  parameters  for  the  lit  region  are  given  by 
2m(Qp)cos0^  ;  p' 
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in  which  p  (Qj^)  is  the  radius  of  curvature  of  the  surface  at  and 


X'  =  2k)!2cos^0^  . 


The  shadow  region  parameters  are  given  by 


r  r  m(t')  .  ,  . 

• 

a  ^ 


u 


ks2  ^ 

~2m(Q^)mfQj^l 


df  ; 


rkpJf) 

;  m(t')  =  — ^2 - 


(A-2d) 


(A-2e;A-2f) 


(A-2g;A-2h) 


It  is  noted  from  Figure  A-1  that  S2  =  (for  0^  =  V2);  82  =  ®o  ~ 

The  functions  F  and  P  are  defined  by  [4J 

.  00  .  2 

F((5)  =  2jjs  f  dx  ;  for  5>0  (for  6<0,  see  Section 

Js  III  of  text).  (A-3) 


(a)  =  ^  7  dx  ; 

fn  -«>  QW2(t) 


1,  for  acoustic  soft 
case 

,  for  acoustic 
hard  case. 


The  a  in  (A-4)  is  positive  in  the  shadow  region;  whereas,  it  is  negative 
in  the  lit  region.  The  shadow  boundary  occurs  at  cfO.  The  Fock  type  Airy 
functions  ^(.j)  and  W^It)  are  defined  in  [P],  such  that 


(b)  SURFACE  DIFFRACTED  RAY  PATH  FOR 
THE  FIELD  POINT  P  IN  THE 
SHADOW  ZONE 


Figure  A-I.  Ray  oaths  associated  with  the  problem  of  scattering  of  an 
obliquely  incident  plane  wave  by  a  smooth  convex  cylinder. 


(A-5a;A-5b) 


?.i  V(t)  =  W^(t)  -  W^fx);  W^rx)  =  —  /  dt 

2 

2 

The  contour  of  integration  r-j  runs  from  ca  to  »,  and  is  the  complex 

conjugate  of  .  The  functions  F(6),  p*(a),  and  q*(0)  are  plotted  in  [4]; 
it  is  noted  that  the  functions  p*(a)  and  q*(a)  are  simply  related  to  P2(a) 

A 

and  P^(<7),  respectively,  via  the  relations  given  in  Equation  (18)  of  Section 
III. 

A  solution  to  the  three-dimensional  (3-D)  problem  of  the  scattering 
of  an  obliquely  incident  plane  wave  by  a  smooth  convex  cylinder  as  in  Figure 
A-1  is  directly  obtained  from  the  solution  in  Equations  (A-la;b)  for  the 
2-D  normal  incidence  plane  wave  case  via  the  method  of  separation  of  variables. 

Let  u^  represent  the  total  exterior  field  for  the  3-D  case,  with  the 
subscript  "3"  being  emploved  to  distinguish  the  3-D  nature  of  the  problem. 

Then,  u.^  satisfies  the  scalar  wave  equation, 

+  k^)  u^  =  0.  (A-F) 


This  obliquely  incident  plane  wave  field  u’  may  be  expressed  in  the  form 


~i  -i 
ul  =  u’^kj  e 


k  =  kcose  I 

r 

k,.  =  ksine 

t  OJ 


where  ^^(k^)  has  the  same  spatial  dependence  as  u^  in  the  2-D  result  of  Eq¬ 
uations  (A-la;b)  except  that  the  "k"  in  those  equations  is  now  replaced  by 
"kj.".  Thus,  U2(k^)  is  the  component  of  the  incident  field  which  propagates 
transverse  to  z  with  a  wavenumber  k^.  One  notes  that  e  is  the  angle  between 
the  incident  ray  direction  and  the  z-axis  as  in  Figure  A-1,  where  the  z-axis 
is  pa>  ellel  to  the  generator  of  the  cylinder.  The  scattered  field  can  also 
be  split  up  as  in  Equation  (A-7),  Consequently,  one  may  write 
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-jk  2 

U3  =  ‘J2(k^)  e 


(A-8) 


where 

(i  *  -  0 


^  3z^ 


(A-9) 


and  ^2  or  ii2(k^)  is  the  same  as  in  Equations  (A-1a;b)  with  the  "k"  in  those 
equations  replaced  by  "k^".  Thus,  is  obtained  via  Equations  {A-la;b), 
(A-7)  and  (A-8)  as 


U3^P)~U3(P)+U3(0p) 


e 


~r 

Pt 


Pt^^sine^ 


-jktil3Sineo-jk^2 

for  P  in  the  lit  region.  (A-lOa) 


U3(P)'-^3(Qa) 


-Jm(Q  )m(Q.  ) 


-jk.t^ 
e  "^sine 


r-jf 


3P-F(x'‘)l+P,(5'’)i 


-jkjSjSinBo-jk^z 

e _  . 

for  P  in  the  shadow  region. 

(A-lOb) 

Ir  obtaining  Equation  (A-lOa)  it  is  assumed  for  the  sake  of  convenience  that 
the  origin  of  the  z  axis  is  located  at  the  point  of  reflection,  Q^;  whereas, 
in  Equation  (A-lOb)  it  is  assumed  that  the  origin  is  at  the  point  of  grazing 
incidence,  Q  .  This  choice  of  the  origin  simply  implies  that  the  axial  (z) 

O 
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separation  between  the  points  Qj^  or  Q^,  and  the  field  point  is  the  distance 
"z".  Thus,  z=J!,2COS0p  in  Equation  (A-lOa)  for  the  lit  zone;  whereas,  z=S2COseQ 
in  Equation  (A-15b)  for  the  shadow  zone.  It  is  easily  verified  that,  the 
parameters  for  the  lit  zone  are 

=  -?mIQj^)f-^/^cos0';  f=l  +  ^  cos^e’;  p[  =  ; 


(A-lla;A-llb;A-llc) 


?.i 


::  2k. (arsine  )  =  2k(i,cos^0\  p"  =  p‘'sin  w  +  p"  cos  u  ; 

Qp  g  T  0  z  0 

(A-lld;A-lle) 


-1-12  -12  ^  ^ 
Ptn  =  cos  u)o  +  sin  o)^;  =  z; 


^  A  A 

Ui  =  U2  X  n. 


(A-llf;A-llg;A-llh) 


The  anqle  of  incidence  {0^)  is  defined  as  usual  by  =  -  s’’'^.n  =  cose\ 

The  angle  (m  )  is  shown  in  Figure  4,  and  p„  and  p  are  the  principal  radii 
O  Z  T 

A  A 

of  curvatures  of  the  surface  along  the  and  directions,  respectively. 

In  the  case  of  a  cylindrical  surface,  p^-w».  The  parameters  p^  and  p^^  are 

A  A  A  A 

the  surface  radii  of  curvatures  in  the  (n,t)  and  (n,bl  planes,  respectively, 

AAA 

where  the  unit  vectors  t,  n,  b  are  shown  in  Figure  4a.  It  is  noted  that 
the  definitions  of  mfQj^)  and  p*"  are  the  same  as  in  Equations  (A-2b)  and  (A- 
2c),  respectively,  except  that  Pg(QR)  in  Equation  (A-lle)  has  the  same  value 
as  pgCQp)  in  Equation  (A-2b)  only  for  8^  =  t\/2  (when  k^=0  and  k^=k);  other¬ 
wise,  Pg(Qp)  of  Equation  {A-2b)  is  really  P^  appearing  in  Equation  (A-lle). 


The  appearing  in  Equation  (A-12a)  is  given  by 


P  =  P^ 


=  (-V)  ^ 

Vsin  0^ 


(A-13) 


The  result  of  Equations  (A-10;b)  or  (A-12a;b)  is  valid  outside  the  paraxial 
regions  of  the  cylinder;  i.e.,  it  is  valid  for  not  close  to  0°  or  180°. 

One  may  view  the  incident  field  u^  as  being  produced  by  a  scalar  point 
source  at  the  point  P  (P=P|^  or  P=Pg)  which  receeds  to  infinity.  If  the  strength 


of  this  source  is  a^,  then: 


“3<'’r>  '  “o  ^ 


-JMsj-sp 

u>)l  -a„c^e  3  3 
lit  zone  ^3 


(A-14a;A-14b) 


/  -j''"3 

“I'V 


•  ''3<'’l|  '=0  ' 

' shadow 


d  ° 


(A-14c;A-14d) 

The  incident  field  uZ  may  be  normalized  to  a  unit  spherical  wave  by  suppressing 
the  factor  c  - ^  for  the  case  of  "plane  wave  incidence". 


Let  a  scalar  point  source  of  strength  b^  placed  at  P  of  Figure  A-1 
generate  an  incident  field  u^  with  the  source  at  P  (i.e.,  P=P|^  or  P=Pg) 
turned  off.  The  total  far  zone  field  of  this  source  in  the  presence  of 
the  cylinder  is  also  obtainable  from  Uj  of  Equations  (A-12a;b)  via  reciprocity 
and  is  given  by 


‘’o  »  i 


P=Pl 

p=p 

^Psl) 


if  P  is  in  the  lit  zone 


if  P  is  in  the  shadow  zone. 

(A-15) 

It  follows  after  some  re-arrangement  of  terms  that: 


#• 


•U 

•J? 


^  T  D-nx'-)]+P  (C*-) 

eJFc*-  I  _ 


p*:  p"" 


■jks- 


^  2  s*^ 

^3 


for  P=P|^  in  the  lit  zone.  (A-16a) 


-  "j^^3  [p 


.  |)-F(X«)].P^{C“), 
^  h 


|dn(Qb) 


-jks^ 


for  P=Pg  in  the  shadow  zone.  (A- 16b) 

In  the  above  solution  of  Equations  fA-16a;b)  for  the  reciprocal  or  the  point 
source  excitation  problem,  all  the  ray  directions  in  Figures  (A-la;b)  must 
be  reversed.  Also,  and  t^  must  be  defined  as 


5“  '  /  2'!^)  df  :  tj-  f  If 


(A-17a;A-17b) 


The  quantities,  p!j^  and  represent  the  reflected  ray  caustic  distances  for 
the  point  source  excitation  case. 


1 


r 


(A-18a;A-18b) 


The  p*^  in  Equation  (A-18a)  is  defined  in  Equation  (A-13).  The  quantities 
p!J^  can  he  obtained  via  Equation  (6a)  of  Section  II-A.  The  quantities  within 
the  radical  signs  outside  the  square  brackets  of  Equation  (A-16b)  are  defined 
as 
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p 

dn(Qu)  s-  /surface  rayr 

w  “  J 


caustic  distance 
p  =  s^+t^  =  for  the  surface  . 

diffracted  ray 


(A-iqa;A-19b) 


It  is  understood  that  the  following  far  zone  approximations  are  implied  in 
Equations  (A-16a;b): 


-jksC  •.  /  i  r.  -jkii.  -jks^ 

^3  ^  ^3 


(A-20a;A-20b;A-20c) 


-^i(Ps)V 


-jks 


3  -jk(s3-S3) 


(A-20d) 


In  conclusion.  Equations  (A-16a;b)  represent  the  far  zone  field  of 
a  scalar  point  source  at  I*”  (see  Figure  la;b)  which  radiates  in  the  presence 
of  a  convex  cylinder.  If  one  suppresses  the  spherical  wave  factor 
c  e  /s^’  in  Equations  (A-16a;b),  one  then  obtains  the  "total  far  field 

pattern"  of  that  point  source  in  the  presence  of  the  cylinder. 
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